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ABSTRACT 

We investigate the stellar populations of 25 massive, galaxies (log[M*/M 0 ] > 10.9) at 1.5 < z < 2 using data 
obtained with the X-band Multi-Object Spectrograph (KMOS) on the ESO VET. Targets were selected to be 
quiescent based on their broadband colors and redshifts using data from the 3D-HST grism survey. The mean 
redshift of our sample is z = 1.75, where KMOS 17-band data probe age- and metallicity-sensitive absorption 
features in the rest-frame optical, including the G band, Ee I, and high-order Balmer lines. Fitting simple 
stellar population models to a stack of our KMOS spectra, we derive a mean age of 1 .03 ;!;q Qg Gyr. We confirm 
previous results suggesting a correlation between color and age for quiescent galaxies, hnding mean ages of 
l.22^lg Gyr and 0.85;!^Qg Gyr for the reddest and bluest galaxies in our sample. Combining our KMOS 
measurements with those obtained from previous studies at 0.2 < z < 2 we hnd evidence for a 2-3 Gyr spread 
in the formation epoch of massive galaxies. At z < 1 the measured stellar ages are consistent with passive 
evolution, while at 1 < z < 2 they appear to saturate at ^1 Gyr, which likely reflects changing demographics of 
the (mean) progenitor population. By comparing to star-formation histories inferred for “normal” star-forming 
galaxies, we show that the timescales required to form massive galaxies at z > 1.5 are consistent with the 
enhanced a-element abundances found in massive local early-type galaxies. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 


1. INTRODUCTION 

At early times, galaxies’ growth is dominated by in-situ 
star-formation fed primarily by direct accretion of gas from 
the cosmic web (e.g., Dekel et al. 2009; Tacconi et al. 2010). 
The relative balance of accretion, star formation, and feed¬ 
back in this “equilibrium growth” phase leads to a tight rela¬ 
tionship between star-formation rate (SFR) and stellar mass 
for normal star-forming galaxies (e.g.. Bower et al. 2006; 
Schaye et al. 2015). The overall normalization of this star¬ 
forming “main sequence” evolves downward from the epoch 
of peak activity at z 2 to the present (e.g., Daddi et al. 2007; 
Rodighiero et al. 2011; Whitaker et al. 2014; Schreiber et al. 
2015), reflecting both a declining accretion rate of gas onto 
halos and the steady depletion of gas via star formation (e.g., 
Lilly et al. 2013). Galaxies appear to spend the majority of 
their lifetimes on the main sequence (e.g., Noeske et al. 2007) 
until their star formation is “quenched” near the Schechter 
mass (e.g., Peng et al. 2010). This quenching process leads to 
formation of a quiescent galaxy population whose subsequent 
evolution is dominated by the assembly of already-formed 
stellar mass via mergers (e.g., Oser et al. 2010; Moster et al. 
2013). 

t Based on observations obtained at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO), Paranal, Chile (ESO program 
IDs 092.A-0091, 093.A-0079, 093.A-0187, and 094.A-0287). This work is 
further based on observations taken by the 3D-HST Treasury Program (GO 
12177 and 12328) with the NASA/ESA Hubble Space Telescope, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 


In the local Universe, quiescent galaxies occupy a well- 
dehned “red sequence” in terms of their color and luminos¬ 
ity (or stellar mass) with a very small intrinsic scatter (cr < 
0.04 mag; e.g.. Bower et al. 1992) suggesting an early, rapid 
formation process which is supported by their old ages and 
high a-element abundance ratios (e.g., Thomas et al. 2005). 
Over the last decade, near-infrared photometric surveys have 
shown that the massive end of the red sequence is in place 
already by at least z ~ 2 (e.g., Williams et al. 2009; Bram- 
mer et al. 2009), consistent with the early formation epochs 
inferred from local elliptical galaxies; however, the details of 
their formation remain poorly understood, motivating a more 
direct study of their stellar populations, and in particular ages 
and a-element abundances, which relate directly to their star- 
formation histories. Only recently-with the development of 
efficient red-sensitive optical and near-infrared detectors-has 
it been possible to measure galaxy stellar populations rou¬ 
tinely up to z 1 (e.g. Sanchez-Blazquez et al. 2009; Belli 
et al. 2015) and as high as z 2 for select sub-samples (Belli 
et al. 2014; Onodera et al. 2014). 

In this Letter we investigate the stellar populations of mas¬ 
sive quiescent galaxies at 1.5 < z < 2 using new data obtained 
with KMOS (Sharpies et al. 2012, 2013) as part of the VLT 
IRIFU Absorption Line GTO survey (VIRIAL; Mendel et al. 
in prep.). These data probe rest-frame V-band absorption fea¬ 
tures for galaxies at z > 1.5 with a resolution of R ~ 3500, 
allowing us to study the stellar content of massive galaxies 
near the epoch when they are being formed. 

Where quoted, magnitudes are on the AB system, and stel- 
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Figure 1. VIRIAL survey sample selection. Light grey points show the dis¬ 
tribution of 3D-HST targets at 1.45 < z < 2 while darker circles highlight 
those galaxies satisfying the selection criteria described in Section 2. Dia¬ 
monds indicate galaxies considered in the present work. The division be¬ 
tween nominally “old” and “young” galaxies based on their colors is adopted 
from Whitaker et al. (2013). The inset histograms show the normalized color 
distribution for the full sample (filled) and the sub-sample used here (open). 
The cun'ent data are roughly representative of the underlying f/VZ-quiescent 
population. 

lar masses have been computed using a Chabrier (2003) stellar 
initial mass function (IMF). We adopt a flat cosmology with 
= 0.7, Hm = 0.3 and Hq = 70 km s“' Mpc“'. 

2. DATA AND ANALYSIS 

VIRIAL survey targets are selected based on combined red- 
shift, magnitude, and color criteria using data from the 3D- 
HST Treasury Program (Brammer et al. 2012; Skelton et al. 
2014) and CANDELS (Grogin et al. 2011; Koekemoer et al. 
2011) in the COSMOS, GOODS-S, and UDS fields. We 
first identify quiescent galaxies in the redshift range 1.45 < 
z < 2 by their rest-frame U-V and V-J colors (hereafter 
UVJ', Figure 1) using the division described by Williams 
et al. (2009, see also Whitaker et al. 2011), supplementing 
publicly-available spectroscopic redshifts with those derived 
using 3D-HST grism and photometric data. We then consider 
only those galaxies with mpMow < 22.5, corresponding to a 
sample which is >90% complete for log(M*/M 0 ) > 10.9. 
The adopted selection includes >95% of the quiescent galax¬ 
ies^ brighter than our adopted magnitude limit. The full 
VIRIAL survey consists of 132 galaxies which satisfy the 
above cuts. In this letter we focus on a sub-sample of the 
spectroscopic data with reliable redshifts obtained during the 
first 2 semesters of observations. These galaxies have a color 
distribution which is roughly representative of the underlying 
target population, as shown by the histograms in Figure 1 . 

2 . 1 . Observations and data reduction 

Observations of the selected VIRIAL targets were obtained 
on 2014 April 6-7, August 18-20, 28-30, September 1, and 

^ Where quiescent galaxies are defined as those with SFRs a factor of 10 
or more below the mean SFR for their stellar mass and redshift given by 
Whitaker et al. (2014). 


November 27-30. Data were taken with a standard object- 
sky-object nodding pattern so that each on-source frame has 
an adjacent sky exposure. We included small pointing offsets 
of 0.1 - 0 . 6 " in each exposure to avoid an undue contribution 
from bad pixels in the final extracted spectra. Typical expo¬ 
sure times were 300s in the YJ band (1-L36p,m), and the total 
integration time per target ranges from 6-11 hours on source. 

Data were reduced using a combination of the Software 
Package for Astronomical Reductions with KMOS pipeline 
tools (SPARK; Davies et al. 2013) and custom Python 
scripts. Dark, arc, and flat frames were processed using stan¬ 
dard SPARK routines. For both object and sky frames, we 
applied a correction for the channel-dependent bias level drift 
determined using reference pixels on the perimeter of each de¬ 
tector prior to reconstructing the data cubes. Cubes were then 
corrected for spatial illumination effects using the observed 
sky-line fluxes on a frame-by-frame basis. Sky subtraction 
was carried out using the methods described by Davies (2007) 
to correct for short-term variability in the OH airglow lines. 
Finally, we corrected for telluric absorption using synthetic at¬ 
mospheric models computed using MOLECFIT (Kausch et al. 
2014). Observations of spectrophotometric standards taken 
on each night were used to rescale the precipitable water va¬ 
por predicted by MOLECFIT in order to better match ob¬ 
served conditions. We then computed telluric corrections for 
each exposure using the predicted time-dependent model at¬ 
mosphere and this derived scaling. 

In order to minimize systematic uncertainties in the inter¬ 
line sky background, one-dimensional spectra for each object 
were extracted using a routine which iteratively fits for an 
exposure- and wavelength-dependent background term con¬ 
temporaneously with extraction of the source spectrum. Since 
our targets are generally undetected in individual KMOS ex¬ 
posures we use the CANDELS/3D-HST F125W image mo¬ 
saics^ (Grogin et al. 2011; Koekemoer et al. 2011; Skelton 
et al. 2014) as a model for source flux distribution, which 
were convolved to match the KMOS PSF determined from 
reference stars assigned to IFUs in each exposure. Uncertain¬ 
ties on the extracted spectra were estimated using 50 bootstrap 
realizations of the input frames. 

2 . 2 . Stellar population measurements 

The signal-to-noise (S/N) of our extracted KMOS spectra is 
generally too low to yield robust stellar population measure¬ 
ments for individual galaxies; however, we can estimate the 
mean properties for galaxies in our sample by stacking their 
spectra given a sufficiently accurate redshift. 

We estimated redshifts by first cross-correlating the ex¬ 
tracted spectra with simple stellar population (SSP) templates 
at a variety of ages. In the current work we limit our sample 
to the 25 galaxies with a clear peak in the cross-correlation 
signal. By requiring well-measured redshifts we impose a 
bias in our sample towards low mass-to-light ratios; how¬ 
ever, the selected galaxies are nevertheless representative of 
the underlying t/y7-selected population (see Figure 1). There 
is generally excellent agreement between our KMOS-derived 
redshifts and those in the 3D-HST target catalog, with a me¬ 
dian absolute deviation of ~1300 km s“\ consistent with the 
uncertainties estimated by Whitaker et al. (2013). Individ¬ 
ual spectra were then normalized using a high-order polyno¬ 
mial before being co-added using inverse-variance weights. 
In total we generated 3 stacks: one for the full sample, and 

^ http :/ /http : //3dhst . research.yale.edu/Data.php 
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Figure 2. Stacked galaxy spectra and their best-fit SSP models for the total (top), red (middle), and blue (bottom) subsamples. Spectra have been normalized by 
a high-order polynomial, and are offset for clarity. Key absorption features are indicated by vertical lines. Shading shows the ±lcr uncertainties on the stacked 
spectra, derived using bootstrap samples. 


two more where we divide the sample into red and blue sub¬ 
populations based on their UVJ colors as in Whitaker et al. 
(2013). The hnal stacks have S/Ns of 16 A“\ 13 A“*, and 
9 and are shown in Figure 2. 

We estimated the properties of our stacked spectra using 
SSP models constructed with v2.4 of the Flexible Stellar Pop¬ 
ulation Synthesis (FSPS) code described by Conroy et al. 
(2009). We ht simultaneously for the SSP-equivalent age t, to¬ 
tal metallicity, and parameters of line-of-sight velocity distri¬ 
bution, using the ensemble sampling code emcee (Foreman- 
Mackey et al. 2013). SSP models were linearly interpolated 
in log t and metallicity on the fly as part of the htting proce¬ 
dure. Fits were limited to A3800^600A in order to exclude 
wavelength regions where only a few galaxies contribute. 

Final measurements of the stellar age were taken as the me¬ 
dian of the marginal posterior distribution; uncertainties were 
estimated using 100 bootstrap realizations of the stacked spec¬ 
tra, recomputing the best-ht age for each using the procedure 
outlined above. We derive a mean age for the full stack of 
1 ■03;!;Qgg Gyr, while for red and blue sub-samples we derive 
ages of 1.22 ;!;q j 9 Gyr and O.SSiliQgj Gyr. In all cases the data 
are consistent with solar metallicity. We stress that these un¬ 
certainties do not include systematic effects due to our choice 
of SSP model. 

3. RESULTS 

3.1. Quiescent galaxy stellar populations af z « 1.75 

The relatively young mean age of massive galaxies in our 
sample suggests that most are recent additions to the quiescent 
population, consistent with other studies that have reported a 
large fraction of post-starburst-like galaxies at high redshift 
(e.g., Whitaker et al. 2012; Bezanson et al. 2013). However, 
the spread in rest-frame colors among quiescent galaxies is 
likely driven at least in part by a variation in their formation 
times or star-formation histories, such that the reddest galax¬ 
ies are also the oldest (Whitaker et al. 2010, 2013). 

In Figure 3 we show the color vs. age derived from our 
KMOS data when they are split into two different bins of 


rest-frame color; the derived ages-as well as the ~0.4 Gyr 
age difference between the two samples-are consistent with 
the evolution of a passively-evolving stellar population (solid 
lines in Figure 3). For comparison, we also consider SSP 
models which host residual star-formation fed by stellar mass 
loss, assuming that 100% of the material lost is recycled into 
stars (dashed lines in Figure 3). While the bluest galaxies are 
compatible with either model, the reddest are unlikely to host 
signihcant ongoing star-formation, more than a few Mq yr“', 
unless they are also preferentially dusty. Of the 25 galaxies 
in our sample, 4 (16%) are detected at 24 p.m in the avail¬ 
able Spitzer-MlPS data, suggesting that they may host ob¬ 
scured star formation; however, there is no clear preference 
for 24 /rm-detected sources to fall in a particular region of the 
t/yy-quiescent selection window (see also Fumagalli et al. 
2014). 

3.2. Mean evolution of red galaxies from z « 2 

In Figure 4 we investigate evolution of galaxy ages from 
z « 2 to the present, incorporating recent results from the lit¬ 
erature. We include age estimates based on spectral htting 
from Choi et al. (2014) at 0.1 < z < 0.7, Belli et al. (2015) 
at 1 < z < 1.5, and Whitaker et al. (2013) at 1.4 < z < 2.2. 
The Belli et al. (2015) data have been converted from their 
T-model parameterization to light-weighted values using the 
age, T, and metallicity as given in their Table 1; we plot here 
the median age and redshift of their data split into two redshift 
bins. Two additional age estimates based on htting to standard 
Lick absorption line indices (e.g. Worthey 1994) from Schi- 
avon et al. (2006) at z = 0.9 and Onodera et al. (2014) at z = 1.6 
are also shown. 

As a comparison to the observed ages we consider the evo¬ 
lution of galaxies which form as a part of the normal star¬ 
forming population, but whose star formation is truncated af¬ 
ter some redshift Zminc so that they evolve passively to the 
present day. We reconstruct mean star-formation histories 
(SFHs) by considering that the evolution of stellar mass is 
driven by a combination of star formation and stellar mass 
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Figure 3. U-V and V-J color as a function of stellar age for the red and 
blue galaxy populations. Lines show predicted color vs. age for SSP models 
at two different metallicities (solid lines), as well for models which include 
star-formation driven by stellar mass loss (dashed lines). Uncertainties show 
the error on the mean derived from bootstrap samples, while thin lines indi¬ 
cate the range of colors present in each sample. Vertical arrows indicate the 
expected change in color for Av = 0.5. The measured ages are consistent with 
expectations from a passively-evolving SED model. 

loss (e.g., Leitner & Kravtsov 2011; Leitner 2012), i.e. 

M*(f) = SFR(M*,z)-3?(f), (1) 

where SFR(M*,z) is the redshift- and mass-dependent SFR, 
and is a convolution of the SFH and stellar mass-loss 
rate, which we estimate using FSPS. Following Lilly et al. 
(2013) we assume that the typical specific SFR (sSFR = SFR 
/ M*) of a star-forming galaxy with mass M* at redshift z is 


sSFR(M*,z) 


0-07(Toiis^)'°'(l+z)' ifz<2 
0-30(p^)'° V+z)^/^ ifz>2. 


( 2 ) 


For a given final stellar mass, M*(ztrunc). we evolve galaxies 
back in time according to Equation 1, using the iterative pro¬ 
cedure described by Leitner & Kravtsov (2011) to solve for 
M*(f), until they reach an initial mass of 10^ Mq; the remain¬ 
ing mass is assumed to form in a burst at z = 8. The resulting 
main-sequence lifetime-from an initial mass of 10*^ Mq to 
4^*(ztrunc)-is Only weakly dependent on the particular value 
of the initial mass, for example adopting 10* Mq results in 
lifetimes which are longer by only ~ 150 Myr. In Figure 4 
we show the evolution of luminosity-weighted age for SFHs 
with 1 < ztrunc < 2.5 and M*(ztrunc) = 1O"M0, as well as for a 


Redshift 



Age of Universe [Gyr] 


Figure 4. Evolution of mean stellar age for massive quiescent galaxies over 
the past ~9 Gyr. Our age estimates, as well as those from Schiavon et al. 
(2006), Onodera et al. (2014), Whitaker et al. (2013), and Choi et al. (2014), 
are measured from stacked data, while the points from Belli et al. (2015) 
are taken as the median of their individually-analyzed galaxy spectra split 
into two redshift bins. We have coiTected the Belli et al. (2015) values from 
their r-model parametrization to light-weighted values for comparison with 
the other data (original values are shown as open diamonds). Error bars in¬ 
dicate either the uncertainty in a given SSP-equivalent age measurement, or 
range in redshift spanned by a given sample. Lines show evolution of the 
light-weighted age for the truncated MS star-formation histories described 
in Section 3.2. The flattening of mean stellar ages at z > 1 likely reflects 
the increased importance of quenching to the formation of massive quiescent 
galaxies at high redshift. 

galaxy which remains star forming until the present day (i.e. 

^trunc “ 0)- 

There are two competing effects which drive the age evo¬ 
lution shown in Figure 4: the quenching of star formation in 
star-forming galaxies, and the passive evolution of already- 
formed stellar mass. At z < 1 the stellar populations of mas¬ 
sive galaxies are consistent with simple passive evolution, 
in agreement with previous spectroscopic studies (e.g., Choi 
et al. 2014), as well as with evolution of the galaxy luminosity 
function over a similar epoch (e.g.. Brown et al. 2007). On the 
other hand, the relatively uniform ages measured at z > 1 sug¬ 
gest that the quiescent galaxy population is being kept young 
(in the mean) by the constant addition of recently-quenched 
galaxies. The transition at z « 1 therefore reflects changing 
demographics of the progenitor population (i.e., progenitor 
bias; van Dokkum & Franx 2001), where the factor of ~3 
increase in quiescent galaxy number density from z = 2 to 
1 (e.g., Tomczak et al. 2014) is dominated by the addition 
of recently-quenched star-forming galaxies, and at z < 1 the 
massive galaxy population grows primarily by the assembly 
of existing quiescent galaxies. Marchesini et al. (2014) ob¬ 
tain qualitatively similar results for the most massive galax¬ 
ies by tracing progenitors at fixed number density (see also 
Papovich et al. 2014). Based on comparison with the mock 
main-sequence SFHs, a ~ 2-3 Gyr spread in formation times 
is required to explain both the relatively old ages found for 
quiescent galaxies already at z = 1.5-2 as well as the ages 
measured at z = 0, but is nonetheless consistent with the small 
scatter observed locally among red-sequence galaxies (e.g.. 
Bower et al. 1992). 
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Figure 5. Stellai* mass growth as a function star-formation timescale for the 
truncated star-formation histories described in Section 3.2. Shaded regions 
trace the mass buildup for galaxies with IO^^Mq < M*(ztrunc) < 10^^'^Mq 
after they join the main sequence population. The point where galaxies have 
formed 50% of their final mass is marked by the dotted line. The vertical 
hatched region shows the range of formation timescales for local early-type 
galaxies implied by their (enhanced) Q-element abundances. The high star- 
formation rates at z > 2 imply that massive quiescent galaxies can form as 
part of the normal star-forming population on very short timescales. 

3.3. The formation timescale of quiescent galaxies 

In the previous Section we considered possible evolutionary 
tracks for quiescent galaxies assuming that their progenitors 
evolve as normal star-forming galaxies. In Figure 5 we show 
the implied mass growth for these models as a function of star- 
formation timescale, in this case defined as the time galaxies 
spend on the main sequence starting from an initial mass of 
lO^M©. 

The chemical abundance patterns of local ellipticals sug¬ 
gest that they form on relatively short timescales, 0.4-2 Gyr 
(e.g., Thomas et al. 2005; Johansson et al. 2012, shown as 
the hatched region in Figure 5), which are consistent with the 
main-sequence lifetimes inferred for massive quiescent galax¬ 
ies at z > 1.5. On the other hand, the rapid evolution of sSFR 
oc (1 -l-z)^ at z < 2 (Equation 2) leads to a dramatic increase 
in the formation time of massive galaxies on the main se¬ 
quence, so that by z « 1 the typical timescale is a factor of 
^2 longer than that inferred from the a-element abundances 
of local early-types. It is important to note that this compar¬ 
ison neglects the role of mergers: Oser et al. (2010) show 
that massive galaxies can accrete ~50% of their final stellar 
mass from z = 1 to 0. Therefore, at least some of the tension 
between main-sequence lifetime and formation timescale im¬ 
plied by Figure 5 for galaxies at z < 1 can be alleviated by 
assuming that the majority of their stars formed in lower-mass 
galaxies, which then assemble at late times. 

4. SUMMARY 

We present new estimates for the ages of massive quies¬ 
cent galaxies derived from deep KMOS observations. These 
data were obtained as part of the on-going VIRIAL GTO 
program, which targets quiescent galaxies selected from the 
3D-HST survey. By stacking the current data we estimate 
a mean luminosity-weighted age of 1.03;!^ Qg Gyr for quies¬ 
cent galaxies at z « 1.75. Separating galaxies based on their 
UVJ colors, we confirm the correlation between color and 
luminosity-weighted age among quiescent galaxies, deriving 


mean ages of 1.22^^^ Gyr and 0.85;|;Qg5 Gyr for the red and 
blue subpopulations. The measured age and color differences 
are consistent with a passively-evolving SSP, suggesting only 
minor contribution from ongoing star-formation. Investigat¬ 
ing the evolution of galaxy ages from z « 2 to 0, we show 
that the current data suggest a spread in formation times of 
at least 2-3 Gyr in order to explain both the mean ages at 
z > 1 and passive evolution at z < 1. There is a clear tran¬ 
sition in the behavior of galaxy ages at z « 1 which likely 
reflects the changing properties of the progenitor population. 
Based on a comparison to the star-formation histories of main 
sequence galaxies at different redshifts, we show that the for¬ 
mation timescale inferred for massive galaxies at z > 1.5 is 
consistent with the enhanced a-element abundances observed 
in local ETGs, suggesting that such systems may plausibly 
form as part of the “normal” star-forming population. 

We are indebted to the entire KMOS instrument and com¬ 
missioning teams for their hard work, which has allowed 
our observing program to be carried out successfully. We 
wish to thank the ESO staff, and in particular the staff at 
Paranal Observatory, for their support during observing runs 
over which the KMOS GTO observations were carried out. 
We also thank the anonymous referee for their constructive 
comments. DJW and ME acknowledge the support of the 
Deutsche Eorschungsgemeinschaft via Project ID 387/1-1. 
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